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ABSTRACT 

We investigate the characteristics of young (< 20 Myr) and bright (Lx > 1 x 10'^^ erg s~^) high- 
mass X-ray binaries (HMXBs) and find the population to be strongly metallicity dependent. We 
separate the model populations among two distinct formation pathways: (1) systems undergoing 
active Roche lobe overflow (RLO), and (2) wind accretion systems with donors in the (super)giant 
(SG) stage, which we find to dominate the HMXB population. We find metallicity to primarily 
affect the number of systems which move through each formation pathway, rather than the observable 
parameters of systems which move through each individual pathway. We discuss the most important 
model parameters affecting the HMXB population at both low and high metallicities. Using these 
results, we show that (1) the population of ultra-luminous X-ray sources can be consistently described 
by very bright HMXBs which undergo stable RLO with mild super-Eddington accretion and (2) the 
HMXB population of the bright starburst galaxy NGC 1569 is likely dominated by one extremely 
metal-poor starburst cluster. 



1. INTRODUCTION AND BACKGROUND 

High-mass X-Ray binaries (HMXBs) form an impor- 
tant observational tool in studies of star-forming regions, 
due to both their relatively high X-Ray luminosity and 
age-specificity, whic h make them good indicators of re- 
cent star formation (jGrimm et al. 2003b). In extragalac- 
tic environments, often only bright HMXBs (here defined 
as Lx > 1 X 10'^^ erg s~^) are detectable. This may sig- 
nificantly affect not only the absolute number of sources, 
but also the time evolution and observable characteristics 
of discoverable HMXBs. 

Recent studies have noted an inverse correlation be- 
tween starburst metallicity and the n umber of observable 
HMXBs produced by each starburst. iMaiid et all (|2004[ ) 
find the Small Magellanic Cloud (SMC) to be overabun- 
dant in HMXBs by a factor of ^^50 compared to the Milky 
Way - though they note that several factors, including (1) 
incompleteness in Milky Way observation, (2) recent star 
formation in the SMC, and (3) the large number of Be- 
HMXBs observed in the SMC, greatly contribute to this 
result. Using population synthesis techn i ques to model 
continuous star formation activity, iDravl (2006) found a 
factor of '--^3 increase in the total number of HMXBs with 
Lx > 1 X 10^^ erg s~^ when the metallicity is decreased 
from solar (Z=Zo) top SMC (Z=0.2 Zq). 

Within this population of bright extragalactic HMXBs 
exists a subpopulation of ultra-luminous X-ray sources 
(ULX; defined here as L^ > 1 x lO^^ erg s-\ the Ed- 
dington limit for a ~ 7 Mq black hole), the exact nature 
of which is not w ell under stood (iFabbiano et al.l llQQTt 
Makishim a et al.l[20 00: Fab biano et al.ll200H iMushotzkvl 
2004t iSwartz et al.l 12004) . Because some of these sys- 
tems have isotropic luminosities above the Eddington lu- 
minosity for black holes (BHs) thought to be produced 
via single star evolution, much work has centered on 
the possibility that ULXs contain an intermediate-mass 



BH gMBH) (IColbert fc Mushotzkvl fTOQl Ivan der Marell 
l2004:lMadhusudhanet al.ll2008[ ). Other arguments point 
to U LXs as a natural ext ension of the HMXB popula- 
tion (jGrimm et al.l 120030 ) . These models may bypass 
the strict isotropic Eddington limit by invokin g photon 
bubb les which can escape from a "leaky disk" (Begelmanl 
120021) . a non-sph erical c omponent of the X-Ray lumi- 
nosity fbeamin g) (lKing , et al. 200 _ ll) or an optically thin 
accre tion disk (lAbramowicz et alj 119881 : lEbisawa et al.l 
I2003D . These stellar mass ULXs may achieve super- 
Eddington luminosities throu gh either strong wi nd accre- 
tion from supergiant donors (jSoria et al.l l2005[) or ther- 
mal timescale Roche-lobe overflow (RLO) (jKing et al.l 
[200?). 

Observations have found a c orrelation between ULXs 
and young starburst activity (jlrwin et al.l 120041 ). The 
most recent observations have suggested that the num- 
ber of ULX sources pe ak between 10-20 M yr after star 
formation has ended (jSwartz et al.l 120091 ). These ob- 
servations have also found ULX to be most prevelant 
in low metallicity environments (f PakuU fc Mirionrl 120021 : 
iZampieri et"all 120041: ISwartz et al...200& ). Several mod- 
els have been posited t o explain the inetalli city depen- 
dence of ULX s ources. 'Ma pelli et al.l (|2009[ ). and later 
Zamp ieri fc Rob erts (2009a) have employed the results of 
evolutionary calculations that indicate an enhanced for- 
mation of high mass black holes (^ 30 Mp,) at low metal- 
licities (jEldridge fc VinkI l2006l : iBelczvnski eFall IMol) . 
The large masses of these BHs provide them with Ed- 
dington limits which reach typical ULX luminosities (^ 
10^° erg s~^). These massive stellar black holes are 
expected to form through di rect stellar collap s e wit n 
no associated natal k ic k (e.g ., iFryer fc Kalogeral (I2001D: 
iMirabel fc Rodriguesl (|2003[ )). However. iZezas et al.l 
(^00__^ found a moderate displacement between ULX 
sources and young stellar clusters in the Antennae galax- 



ies, suggesting that the ULX population received a natal 
kick, indicating that these systems either did not un- 
dergo direct collapse, or that direct collapse events are 
still associated with a small to moderate natal kick. 

Seeking to bridge the gap between theoretical predic- 
tions and observation, in this paper we analyze the for- 
mation of bright HMXBs in the first 20 Myr after star- 
burst activity at several different metallicities. We show 
that, while our models can produce bright HMXBs at 
all metallicities, the number of systems, time evolution, 
and orbital period distribution are strongly metallicity 
dependent. In particular, while our research supports 
the observation that ULX are preferentially formed in 
low-metallicity environments, we find that this relation 
is not related to the formation of large (~ 30 Mq) direct- 
collapse BHs at low metallicities, but instead to the dy- 
namics of the pathway which allows RLO onto black 
holes of ~10-15 M0 BHs, with mild super-Eddington lu- 
minosities. 

In Section[2]we describe our simulation codes and mod- 
eling assumptions. In Sections [3] and |4] we analyze our 
results on the properties of bright HMXBs and discuss 
their physical origin in the context of their evolutionary 
formation pathways. In Section [5] we discuss the rela- 
tions between our results and previous work, specifically 
noting the results for the ULX population in Section [5H 
Further applying our results to observed systems, in Sec- 
tion l5.2l we discuss the production of HMXBs in the star- 
burst galaxy NGC 1569 in the context of the ages and 
metallicities of the known super star clusters. Finally, 
we summarize our conclusions in Section [Sj 

2. SIMULATION CODE AND MODELS 

Employing a sophisticated and recently updated popu- 
lation synthesis code, StarTrack (Bclczynski et al. 2002, 
120081 ). we generate HMXBs from a variety of input pa- 
rameters using Monte Carlo methods. We note that the 
parameter space for Monte Carlo population synthesis 
of HMXB systems is formidable in size: thus a full ex- 
ploration of the parameter space is not computational 
feasible. We take default valu es of many parameters, fol- 
lowing the layout set forth in iBelczvnski et al.l ()2008[ ) . 

For our initial conditions, we employ a Salpeter 
(M~^'^^) initial mass function, setting the minimum pri- 
mary mass at 4 Mq and choosing the secondary mass 
via a flat mass ratio distribution with a minimum sec- 
ondary ma ss of 0.08 Mq. Following the prescription of 
lAbtl ()198tl ). we create a distribution of initial binary sep- 
arations which is flat in the logarithm with an upper limit 
of 10^ Rq and a lower limit such that the primary star 
initially fills at most half of its Roche Lobe. We employ 
a thermal distribution for initial eccentricities (Heggie 

[19751). 

The StarTrack code then models the evolution of 
each individual star using the stellar tracks modified 
from those used in the Single Star Ev olution (SSE) 
code (|Pols et all 1 19981: iHurlev et al][2000l ). These mod- 
els include metallicity dependent evolution parameters 
which affect the radius, temperature, and wind mass loss 
rates of each star in our simulation. We also employ 
a revised model for win d accretion from massive stars 
(jBelczvnski et al.l ()2010r), which employs m odifications to 
adapt for the resuhs of lVink et air()2001[ )). 

The stellar characteristics of each star are modified 



through their binary interactions. The most important 
stellar effects come from stable mass transfer through 
Roche Lobe Overflow (RLO) and unstable common enve- 
lope (CE) phases. Stable RLO phases can be treated by 
comparing the corresponding timescale for mass transfer 
from angular momentum losses and nuclear evolution of 
the donor with the donor's thermal timescale. As long as 
the donor remains in the rmal equilibrium the ma ss trans- 
fer rate is then stable (jBelczynski et al.l 120081 Section 
5.1). C E phases are treated using the energy formalism 
(jWebb ink 1984), where the binary is immediately circu- 
larized and the amount of energy needed to remove the 
envelope of the primary star is subtracted from the po- 
tential energy of the binary orbit. If this process results 
in a stable binary system, we continue to track the evolu- 
tion of the binary system. However, if there is not enough 
orbital energy in the binary to remove the envelope with- 
out creating a merger of the two cores, we immediately 
end the simulation. For main sequence (MS) donors, the 
common envelope is invariably assumed to lead to a bi- 
nary merger. Hertzsprung gap (HG) stars have not yet 
fully formed clear core-envelope structure and the inspi- 
ral process in the CE phase is ve ry likely to lead to a bi- 
nary merger fjTaam fc Sandquistl f2000). However, due to 
the uncertainties in the common envelope input physics, 
we either assume that all binaries merge in CE initiated 
by HG star (our default models) or are allowed to sur- 
vive provided that there is enough orbital energy to eject 
envelope (alternative model). For further details on the 
importance of this parameter choice to the formation of 
binary systems, see IBelczvnski et all ()2007|) . 

We set the maximum neutron star (NS) mass at 
2.5 Mq and the minimu m CO core mass for direct col- 
lapse BHs to be 7.6 Mq (jFrveJ 119991 : iFrver fc Kalogeral 
l2001h . At the time of compact-object formation, we as- 
sign a random natal kick drawn from a single Maxwellian 
distribution with a dispersion velocity of 265 km s~^ 
(jHobbs et al.l 120051 ). We scale down the natal kick as- 
signed to BHs by multiplying the kick by the fraction 
of material which does not fall back onto the compact 
object. In situations where BHs form through direct col- 
lapse (e.g., where the entirety of the supernova ejecta 
falls back onto the BH and there is no mass loss other 
than the BH binding energy), we apply no natal kick in 
our default models. 

We note that a ULX cutoff of Lx = 1 x 10^^ erg 
s^^ sets an Eddington-limited BH mass of ~7 Mq. 
However, several models allow ing super-Eddington ac- 
cretion have been postulated (lAbramowicz et al.l 119881 : 
iKing et all 120011 : iBegelmanI I2002D . Instead of assuming 
a specific physical model for super-Eddington emission, 
we simply set the luminosity to be the lesser of (1) the 
luminosity calculated without regard to the Eddington 
limit or (2) ten times the Eddington limited luminosity 
for BH accretors and two times the Eddington limited 
luminosity for NS accretors. Such super-Eddington lumi- 
nosities for BH s have been found in simulations of slim- 
accre tion disks fAbramowicz et al.'1988:'Begel manll2002l : 
Ebisawa et al.l 12003; So crates fc Davis 2006) and pho- 
ton bubble instability (jRuszkowski fc Begelmanl |2003| ) 
These effects may even allow violation s of the Edding- 
ton l imit by a factor as large as 100 (jBegelmanl 120021 
l2006f ) so our super-Eddington model may remain conser- 
vative. On the other hand, NS accretors seem to remain 



within twice the Edding ton rate (jLevine et aLl I1991L 
1 19931 [Grimm et aT]|2003aD . This approximation has sub- 
sequently been employed in several theoret ical models 
of luminous X-Ray bina ry formation (Rap paport et aLl 
l2005t iFragos et alJl2008l ) , and we adopt this assumption 
throughout. In order to compare our results with Chan- 
dra observations, we apply an estimated correction for 
the Chandra energy band as describ ed and justified in 
Section 9.1 of iBelczvnski et all ()2008[ ). 

We make two important parameter space choices 
throughout this work. We first choose to evaluate only 
systems younger than 20 Myr of age. This cutoff corre- 
sponds both to the first appeara nces of electron-c apture 
supernovae (ECS) (|Nomotoll98llPoelarends et al.ll2008D 
in our simulations as well a s the last iron-core collapse 
SN. In iLinden et all (|2009( ). we have found that ECS 
are highly associated with the production of Be-HMXBs. 
Reg ardless of the accuracy o f this association, we note 
that iMcSwain fc GiesI ()20Q5[ ) found the observed num- 
ber of Be-HMXBs to be greatest between 25-100 Myr 
after star formation, relegating the vast majority of Be- 
HMXB activity until after the end of the simulations and 
the questions considered here. Be-HMXBs have been 
found to dominate well studied systems, such as the 
SMC, with transi ent luminosities exceeding our brigh t 
luminosity cutoff ()Liu et al.ll2000HAntoniou et al.ll2008[ ). 
However, current models of non-spherical winds from Be 
donors are not well refined, and thus the number and du- 
ration of bright Bc-HMXB activity is a major source of 
uncertainty for our models. Secondly, we have chosen a 
luminosity cutoff of Lx = 1 x 10^^ erg s~^, as this is the 
sensitivity limit of many Chandra observations of young 
starbursts and HMXB populations. 

Another possible mechanism for the production of 
bright HMXBs are short-live d phases of unstable atmo- 
spheric Roche Lobe overfiow ()Savoniielll979l Il983| ) . Due 
to the very short timescale of atmospheric RLO (on the 
order of 10^ yr (iSavoniid 119831) compared to 1 0^ yr for 
thermal timescale RLO (|Belczvnski et al.l 120081 ) ). we do 
not consider these systems in this work. However, in or- 
der to estimate the relative error caused by disregarding 
these systems, we investigate the percentage of atmo- 
spheric RLO events which result in stable RLO. Given a 
Salpeter initial mass function, the average initial mass of 
SN progenitors in our simulation is ^24 M©. BHs formed 
form these primary stars will lose approximately 2/3 of 
their initial mass before CO creation. Since we employ a 
flat secondary mass distribution, we then expect at least 
1/3 of initial secondaries to be less massive than the CO 
- and thus to undergo a thermal timescale RLO-HMXB 
period after atmospheric RLO. Thus, we expect the omis- 
sion of atmospheric RLO to create at most a 20% under- 
estimate in the population of bright RLO-HMXBs in our 
simulations. Since most RLO-HMXBs already accrete at 
or above the Eddington limit, we expect no increase in 
ULX activity from atmospheric RLO. 

To asses the major uncertainties of single and binary 
evolution, we evaluate changes due a number of impor- 
tant parameters individually. The parameters most likely 
to effect HMXB production include: CE efficiency, the 
treatment of HG CEs, natal kick distributions, the exis- 
tence of natal kicks for direct collapse SN, and violations 
of the Eddington limit. The exact role of each of these 
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(ULX, bottom) over the first 20 Myr after a starburst of 10'' 
Mq. Numbers are generated as the average of ~122 simulations. 
(Default evolutionary model, lOx Eddington limit hx allowed) 

parameters will be examined in Section |4l 

3. RESULTS FOR DEFAULT MODEL 

Our models show conclusively that the amplitude 
(number of systems) and time dependence of HMXB ac- 
tivity are strongly metallicity dependent. In Figure (TJ we 
show the number of HMXBs with Lx > 1 x 10^^ erg s~^ 
(bright, top) and Lx > 1 x 10^^ erg s^^ (ULX, bot- 
tom) as a function of time for a starburst of 10^ M© at 
metallicities of Z=(Zq, 0.4Zq, 0.2Zo, O.OSZq, O.O2Z0)0 
All computational simulations presented throughout this 
paper were done for an initial population of 10^ binaries 
with primary mass greater than 4M0, which making as- 
sumptions for a binary fraction of 1.0 and a Salpeter 
IMF, translates to approximately 122 realizations. In all 
plots we have shown the average of these realizations. 
Evaluating our population of bright sources over the en- 
tire 20 Myr simulation, we find that decreasing metallic- 
ity monotonically increases the number of HMXBs by 
a factor of ^3.5 between 7j—7jq and Z=0.02 Z0. At 
the ULX luminosity cutoff (bottom), the low metallic- 
ity sources are even more dominant, with a factor of '^5 
change from our lowest to highest metallicity and with a 
increase of ~2.5 between Z=Zq and Z=0.4 Z©. 

The time evolution of our HMXB population is also de- 
pendent on the starburst metallicity. During the period 
between 5-10 Myr after star formation, the metallicity 
dependence of bright HMXB is non-monotonic, with a 
peak at Z=0.4 Z© and less than half as many systems at 
both metallicity extremes. The reasons behind this non- 
monotonicity are discussed in Section 14.21 This trend 
disappears after « 10 Myr, after which low metallicity 

^ Throughout this paper, we generally refer to systems of Z0 
and O.4Z0 as high metallicity environments, systems with O.2Z0 
as a medium metallicity environment, and systems with O.O5Z0 
and O.O2Z0 as low metallicity environments. This naming scheme 
is based on the trends we see in the HMXB population as a function 
of metallicity, but may be at odds with observational conventions. 



oa 

X 

2 
X 

x: 

BO 




0.1 



0.1 



M iiii| 1 — I I I Mill 1 — I I I mil 1 — I I I mil 1 — I I 1 1 1 li 

Z=0.4 Zg 

ivjsv Z=0.2 Zg 

\| — -Z=0.05 Zg 

- Z=0.02 Z, 




I I ' I ' I " = "■ I 



1036 1037 1038 1039 IQW 

X-Ray Luminosity Function (erg s"' 



lO" 



Fig. 2. — Cumulative X-Ray Luminosity functions for HMXBs 
at Z=(Z0, O.4Z0, 0.2Zq, O.OSZ©, O.O2Z0) wlien strictly applying 
the Eddington limit (top) and for models allowing super-Eddington 
accretion using our default formalism (bottom). Results are taken 
over the entire 20 Myr simulation and are normalized to the total 
number of systems with hx > 1x10^® erg s~^ at each individual 
metallicity when super-Eddington accretion is enabled. (Default 
evolutionary model) 



systems dominate the bright HMXB population by a fac- 
tor of 5. We note that while the number of low metal- 
licity systems declines after 20 Myr, they continue to 
dominate over the high metallicity population, creating 
an even stronger trend towards low metallicity HMXBs 
when we continue our simulation to later time periods. 

At the ULX cutoff (bottom), we note that after 10 Myr, 
the high metallicity HMXB population almost entirely 
disappears, while the low metallicity populations are 
reaching their peak. In this period, we find HMXBs at 
Z=0.02 Z0 to outnumber Z© HMXBs by a ratio of nearly 
1000 to I (with considerable statistical uncertainties due 
to the very small number of high metallicity systems pro- 
duced). This lies in stark contrast to the approximately 
2 to 1 ratio found at our bright luminosity cutoff. 

We note that our determination of HMXB numbers 
could be greatly influenced by our choice to allow mild 
violations of the Eddington limit. In Figure [2] we show 
the X-Ray luminosity function (XLF) at each metallicity, 
given both strict adherence to the Eddington cutoff (top) 
as well as our super-Eddington formalism (bottom). We 
normalize the number of HMXB at each metallicity to 
the number above our bright cutoff in our default model. 
We note that the XLFs in our study are consistent with 
power-laws up to a luminosity of hx = 1 x 10'^^ erg s~^ 
and are not metallicity dependent, consistent with ob- 
servations of local galaxies. We caution, however, that 
our modeled XLF should not be considered as a model 
for the global HMXB population in galaxies with young 
or continuous star formation, for several reasons: (1) we 
consider only HMXBs formed within the first 20 Myr 
after star formation, while all observable populations 
contain both young and medium age HMXBs, (2) we 
plot only the orbit averaged luminosity due to spheri- 
cal winds, while the majority of observed extragalactic 
systems are found during burst activity, and (3) obser- 
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Fig. 3. — Cumulative distribution of compact object masses for 
sources at Z=(Z0, O.4Z0, O.2Z0, O.O5Z0, O.O2Z0) for systems with 
Lx > lO^^ergs-i (top) and Lx > lO^^ergs"^ (bottom). All re- 
sults are normalized to the total number of systems throughout 
each simulation, which here corresponds to 122 realizations of lO'' 
M0 starbursts with an average HMXB number (over 20 Myr) of 
5.73 at Z=O.O2Z0, 5.25 at Z=O.O5Z0, 3.81 at Z=O.2Z0, 3.79 at 
Z=O.O4Z0, and 1.71 at Z=Z0. The number of systems with com- 
pact object masses greater than 40 M0 is less than 5 % at all 
metallicity and luminosity cuts. We note that the number of sys- 
tems at ULX cutoff and Z=Z0, is very small, and this result is 
statistically uncertain. (Default evolutionary model, lOx Edding- 
ton hx allowed) 



vatio ns of the extragalactic XLF (see e.g. IGilfanov et al.l 
l2004r ) are highly dependent on the assumed SFR in each 
galaxy, and often have statistical uncertainties due to 
small number statistics. Given these caveats, careful 
analysis of the XLF shape and of comparisons to ob- 
servational XLFs are beyond the scope of this present 
study, and thus we do not expand further. Such an anal- 
ysis would require developing models which target spe- 
cific known galaxies, using their star formation history 
as model input. Here we instead focus on the bright- 
est HMXBs in starburst environments and how metal- 
licity influences their formation pathways. Finally, we 
note that the simple power-law explanation underscores 
the fact that our resultant population should not be 
highly sensitive to our choice of a luminosity cutoff at 
Lx = I X 10'^^ erg s~^ - a result which will become clearer 
as we analyze the pathways through which HMXBs form. 
At all metallicities, we find that our choice to allow 
super-Eddington accretion increases the total number 
of systems above our bright luminosity cutoff by less 
than 2%. However, Figure [5] also shows that allowing 
super-Eddington accretion greatly affects the number of 
systems above our ULX cutoff. In the Eddington lim- 
ited case (top), we note a sharp cutoff in the XLF at 
L:r > 2 X 1039 erg s-\ with nearly 30% of HMXB having 
luminosities 1 x 10^^ erg s""'^ < L^, < 2 x lO'^^ erg s~^. The 
proximity of these luminosities to our ULX cutoff is not 
surprising, as we have defined the ULX cutoff to be near 
the Eddington luminosity for H accretion from a 7 M© 
BH. However we not e that no such cutoff is observed by 
IGilf anov et al.l (12004'), which is at odds with simulations 
finding such a large fraction of bright HMXBs radiating 
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lO'^^ergs"^. (Default evolutionary model, lOx Eddington hx al- 
lowed) 



at the Eddington limit. When we ahow super-Eddington 
accretion (bottom), we note that between 20-30% of our 
bright HMXB population now has a luminosity greater 
than the L^; > 2 x 10^^ erg s~^ cutoff, with a new power 
law of approximately 0.75. This implies that many of 
our Eddington-limitcd systems have mass transfer rates 
capable of super-Eddington accretion if the accretors are 
able to accept this additional mass. 

In order to further analyze the observable parame- 
ter space of our bright HMXBs, in Figure |3] we plot 
the masses of the compact objects (COs) in our mod- 
els, noting several important results. First, we see very 
few NSs (systems of 1-2 M0) at our bright luminosity 
cutoff (< 5%) and no NSs our our ULX cutoff. More 
importantly, the ULX population does not have signif- 
icantly more massive accretors than the bright HMXB 
population. At low metallicity, the vast majority of 
ULX accretors lie between 10-15 M©. We note that 
although simulations show that low metallicities tend 
to favor la rge single-star BH masses (jBelczynski et al.) 
120061 120091 ) , the BH masses we observe in low metallic- 
ity environments does not approach the maximum possi- 
ble masses from single star models. Instead we find that 
the largest BH masses correspond to moderate metallic- 
ities (Z=0.2 Z0). This absence of massive ULX accre- 
tors lies in direct contrast to earlier discussions of reasons 
for incr eased ULX production in low-metallicity environ - 
ments ()Mapelh et al.l 120091: [Zampieri fc Roberts! I2009bl ). 
We will further analyze this issue in Section [STTl 

These results demonstrate the substantial variation be- 
tween bright and ULX HMXBs, as well as high and low 
metallicity binaries. The metallicity dependent time evo- 
lution at both luminosity cuts in our simulation indicates 
that we are likely dealing with an HMXB population 
which is composed of several distinct sub-groups. In or- 
der to analyze the trends which affect our bright HMXB 
population, we will now examine the evolution pathways 
which are responsible for both bright and ULX HMXB 
activity. 



4. PATHWAYS AND MODEL VARIATION 

We noted several interesting trends in Section |3] which 
indicate that our population of HMXBs may not have 
uniform characteristics. However, the smoking-gun 
demonstrating the metallicity dependence of the bright 
HMXB population is shown in Figure |4l where we plot 
the cumulative period distribution of HMXBs at each 
of our test metallicities. We note that there is a stark 
contrast between the orbital periods of the HMXB pop- 
ulation at Z=0.02 Zq and Z=0.05 Z0 and the popula- 
tion at Z—OA Zq and Z=Zq, with Z=0.2 Z0 acting as 
an intermediate case. For low metallicity systems, the 
vast majority (>60%) have periods of less than 3 days, 
while 90% of high metallicity systems have a period of 
greater than 10,000 days. At the Z=0.2 Z0 metallicity, 
approximately 20% of the systems have periods of less 
than 1 day, with the remaining 70% of systems having 
periods of more than 10,000 days. We note that at high 
metallicities, there are almost no intermediate systems of 
periods between 1-10,000 days, while at low metallicity, 
this accounts for approximately 20% of bright HMXB. 

In order to isolate the effect of metallicity on our 
HMXB population, we analyze our model populations 
of bright HMXBs and identify six distinct binary evo- 
lution pathways: (1) systems undergoing stable mass 
transfer via Roche lobe overflow (RLO-HMXB), (2) sys- 
tems with wind accretion from a MS donor that has 
moved through a CE phase, (3) systems with wind ac- 
cretion from a MS donor that has not moved through a 
CE phase, (4) systems with wind accretion from a (su- 
per)giant donor (SG-HMXB), (5) systems with wind ac- 
cretion from a He-rich donor that has moved through a 
CE phase, and (6) systems with wind accretion from a 
He-rich donor that has not moved through a CE phase. 
Given these conditions, each bright HMXB falls into ex- 
actly one of these catagories at each moment in time, al- 
though HMXBs can move between different populations 
over time. We note that thes e pathways are very simi- 
lar to those flrst described by iVan Bever fc VanbeverenI 
(|2000f) in their study of the important HMXB contribu- 
tions to the galactic X-Ray luminosity function (XLF), 
with agreement on the dominant pathways in the high lu- 
minosity regime. Throughout the rest of this section, we 
will restrict ourselves to examining metallicities of Z^Z©, 
Z=0.2 Zq and Z=0.02 Z© for simplicity, as we find these 
metallicities to adequately demonstrate the main effects 
and dependencies of our models. 

Due to the high luminosity cutoff of 1 x 10"^® erg s^^, 
our simulations do not find any wind accretion systems 
with main sequence or Helium rich donors. We note that 
the spherical winds of MS and He-rich stars are not gen- 
erally strong enough to promote such extreme HMXB 
luminosities. However non-spherical winds (e.g., due to 
rotation) may create bright systems which our models 
would not account for. Specifically, observations point 
towards the existence of a large population of tra nsient 
HMXBs with Be donors ()Liu et al. 2000 : Antoniou et al.l 
120081 ) that can temporarily exceed our luminosity cut- 
off. However, as noted in Section [21 we have chosen here 
to study only systems younger than 20 Myr, where we 
believe the observed sample to contain few Be-HMXBs. 

Thus, we find our bright HMXB population to be 
dominated by systems moving through either the RLO- 
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Fig. 5. — Categorization of HMXBs which shows the (Column 1) time dependent HMXB number (2) cumulative X-Ray luminosity function 
(3) cumulative orbital period distribution (days), (4) cumulative eccentricity distribution (5) cumulative initial mass ratio distribution 
(secondary mass/primary mass), (6) cumulative initial periastron separation (in R©) for systems from (Row 1) RLO-HMXBs (2) SG- 
HMXBs. Each plot except for the first is normalized to the number of systems moving through each pathway at each given metallicity. For 
readability, data are shown only at reference metallicities of Z=Zq, Z=0.2 Zq, Z=0.02 Zq. (Default evolutionary model, lOx Eddington 
hx allowed) 

HMXB or the SG-HMXB pathways. In order to bet- 
ter understand the intial and final conditions which cre- 
ate systems of each class, in Figure [S] we show vari- 
ous observable parameters and initial conditions for sys- 
tems moving through each of our two dominant path- 
ways. We immediately note two important results: 
(1) low-metallicity HMXBs are dominated by systems 
moving through the RLO-HMXB pathway, while high- 
metallicity HMXBs are dominated by systems moving 
through the SG-HMXB pathway; (2) to a rough approx- 
imation, metallicity does not greatly affect the observable 
properties of systems moving through a given pathway, 
instead it affects the number of HMXBs in that pathway. 
In order to better understand the physics which controls 
each HMXB pathway, we will individually examine each 
channel and explain the physical properties that create 
the observable HMXB population from the typical pro- 
genitor parameter space for each channel. 
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4.1. RLO-HMXB Pathway 

Roche lobe overflow provides an important method for 
creating bright HMXBs (with typical luminosities above 
Lx > 1 X lO^^* erg s^^) due to its effectiveness in trans- 
ferring matter onto the compact object accretor. In or- 
der to obtain stable mass transfer through RLO, these 
HMXBs must have both short orbital periods (often < 
1 day. Figure [5l row 1, column 3), and a mass ratio that 
is near unity during the HMXB phase. Due to the exis- 
tence of stable mass transfer, these systems are assumed 
to be circulariz ed in our models (for a discussion of this 
assumption, see lSepinskv et all (|2007[ )V 

In order to create HMXBs with small orbital periods, 
all RLO-HMXBs evolve through a CE phase prior to 
the supernova of the primary star. CE phases are nec- 
essary for RLO-HMXB progenitors because the orbital 
separation necessary to survive a strong SN natal kick 
(~ 10 Kq) is substantially less than the orbital separa- 
tion necessary to prevent the binary merger of unevolved 
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Fig. 6. — Average number of RLO-HMXBs with Lx > 
10'^^ ergs"-"^ over the first 20 Myr after a starburst of 10^ Mq. for 
Maxwellian natal kick distributions with different dispersion veloc- 
ities (Natal kick velocity varies from default model, lOx Eddington 
hx allowed) 



massive stars (> 100 R© at all metallicities). Further- 
more, a very tight binary system (typically < 5 R©) is 
necessary to allow eventual RLO of the donor star. Be- 
cause the creation of CEs requires a large mass ratio be- 
tween the primary and secondary stars (Figure [5l row 1, 
column 5) the majority of RLO-HMXB systems (> 90% 
at all metallicities) contain an unevolved donor accreting 
onto the compact object. 

The dominance of the CE driven pathway explains the 
time evolution of RLO-HMXBs. We do not begin to see 
bright systems until after 6 Myr, the lifetime of the heav- 
iest stars which are able to enter the (super) giant phase 
(and thus undergo survivable CE evolution according to 
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Fig. 7. — Maximum radius of stars during the HG (solid black) 
and (super)giant (dashed red) st ages for stars of vary ing metallicity, 
as determined by the SSE code l IHurlev et al.|[2000f) 



(jTaam fc Sandauistl[2000l) ). The number of bright RLO- 
HMXBs continues to increase as the number of primary 
supernovae increase, since no evolution of the donor star 
is required. This upward trend ends at « 13 Myr, corre- 
sponding to the latest forming supernovae in our simula- 
tion. The number of RLO-HMXBs exponentially decays 
after 20 Myr as the donors begin to move off of the Main 
Sequence. This characteristic evolution is observed at all 
metallicities. 

However, the lack of donor evolution introduces a prob- 
lem for the creation of RLO-HMXBs. The end product 
of CE evolution should be a He-rich primary alongside a 
MS secondary. The He-rich primary is always at least as 
massive as the eventual CO. However, the system must 
not be in RLO at the end of the CE phase, or the sys- 
tem would be considered a merger via the energy for- 
malism. Since the evolution of the donor star is mini- 
mal, another mechanism must be responsible for moving 
the He-rich/MS system into a CO/MS system in RLO. 
This extra force is provided by the natal kick provided 
to the system at CO formation. While the natal kick 
will likely add energy to the binary system, it's more 
important contribution is to take the circularized post- 
CE binary, and move it into a highly eccentric orbit. 
Even a more distantly separated system is thus likely to 
be in RLO during periastron - and the frictional forces 
of RLO are then assumed to circularize the orbit and 
create a tightly bou nd system in perpetual RLO (see 
ISepinsky et al.l [20071 ). In Figure El we show the num- 
ber of RLO-HMXBs for different Maxwellian dispersion 
velocities for the natal kicks imparted to our CO pop- 
ulation. We note that the number of systems increases 
significantly when larger natal kicks are employed. At 
high metallicities, almost no RLO-HMXBs exist when 
natal kicks are disabled. We note a similar correlation 
between natal kicks and the number of RLO X-Ray Bi- 
naries has been previ ously theorized in the case of low 
mass X-Ray binaries (jKalogera fc Webbinklll998l ). 

At our Z=0.02 Zq metallicity, we note a persistent class 
of HMXBs which exist even when natal kicks are dis- 
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Fig. 8.— Average number of HMXBs at Z=(Z0, O.2Z0, O.O2Z0) 
for Lx > 1 X 10'^^ erg s~^ (bright, top) and Lx > 1 x 10'^^ erg s"^ 
(ULX, bottom) over the first 20 Myr after a starburst of 10^ Mq for 
systems going through RLO-HMXB phases with Hertzsprung Gap 
common envelope evolution enabled using the energy formalism. 
(HG CEs allowed, lOx Eddington Lx allowed) 

abled. This residual, which contains almost 50% of the 
total number of systems, is composed of compact objects 
resulting from the direct collapse of particularly massive 
BH progenitors into BHs of 10-15 Mq. Interestingly, we 
find this cl ass of systems to ex ist only when the wind pre- 
scription of lVink et al.l (120011) is emp loyed. Compared to 
the formulas of (jHurlev et al1l200Cl[) . the Vink formula- 
tion predicts significantly lower wind mass loss in the low 
metallicity regime. This allows significantly less massive 
ZAMS stars to end their lives as fallback BHs, compared 
to systems in the high metallicity regime. These ZAMS 
stars can be small enough to move through a SG phase 
and form a tight binary orbit with a MS donor. This 
region of parameter space is not available when either 
stronger wind models are used or when the stellar metal- 
licity is high, as the minimum ZAMS mass necessary 
to create a fallback BH is then larger than the maxi- 
mum ZAMS mass necessary to move through a stable 
(super)giant phase. Upon the formation of the BH, the 
small mass loss from the binding energy of the progenitor 
star creates an eccentricity that brings the system into 
RLOF and creates a bright HMXBs. 

The dominance of the CE based pathway for RLO- 
HMXB formation explains the inverse correlation be- 
tween RLO-HMXB number and metallicity. In order 
for a CE to form, the Roche Lobe of the primary star 
must lie in between the maximum stellar radius during 
the HG and the maximum stellar radius during the SG 
phase. If the Roche lobe is smaller than the maximum 
HG radius, the system will undergo a CE during the HG 
phase which we assume leads to binary merger. If the 
Roche lobe is larger than the maximum SG radius, then 
the system will never enter a CE phase, and thus will not 
form a RLO-HMXB. The survival of high mass binaries 
through the common env elope phase at low met allicity 
was explained in detail bv iBelczvnski et all (|2010[ ) in the 
context of double compact object formation. 

In Figure [71 we show an illustrative plot from the Sin- 
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Fig. 9. — Average number of RLO-HMXBs with Lx > 
10^''ergs~l over the first 20 Myr after a starburst of lO'' M©. 
for different values of the Common Envelope efficiency 



gle Star Evolution code (jHurlev et al.ll2000[ ) of both the 
maximum HG radius (ktype = 2 in SSE) and the maxi- 
mum SG radius (ktype = 4) as a function of metalhcity 
for star with a Zero Age Main Sequence mass of 30 M0 , 
which is typical for the primary star in RLO-HMXB pro- 
genitors. We note that while the maximum SG radius 
remains large even at low metallicities, the maximum 
radius during the HG phase drops precipitously with de- 
creasing metalhcity. These curves place strict upper and 
lower limits on the periastron Roche lobe separation of 
the binary system. We note that our initial parameter 
space of orbital se parations is set to be flat in the log- 
arithm ()Abtl 119831) . Since the necessary binary separa- 
tion is simply a multiplicative constant times the Roche 
Lobe radius (depending on the mass ratio of the system), 
the parameter space of allowable stable SG CEs is much 
wider at low metalhcity than at near solar metalhcity. 
While not all of these CEs are survivable via the energy 
formalism, the trend still creates a large preference to- 
wards low metalhcity RLO-HMXBs. 

Since HG CEs eliminate a large fraction of the parame- 
ter space for RLO-HMXB formation at high metallicities, 
it is worthwhile to determine the effect of removing this 
restriction. In Figure El we investigate the importance of 
this assumption by showing the number of bright (top) 
and ULX (bottom) RLO-HMXBs as a function of time 
for a simulation where we allow CEs in the HG phase to 
survive if they are energetically stable. We note two re- 
sults: (1) the number of RLO-HMXB systems increases 
at all metallicities, as all systems with a HG primary 
were previously assumed to merge, while the new analy- 
sis instead allows some to survive and (2) the first RLO- 
HMXB systems form at much earlier time periods than 
in our default simulation, as primaries too massive to 
move through a (super)giant phase can now go through 
a stable common envelope phase before the 6 Myr genesis 
of (super)giant systems. 

The effects of allowing HG CEs can be drastic, espe- 
cially at high metallicities, where so few RLO-HMXB 
progenitors survive the CE during the SG phase. At 



Z=Z0, we see a factor of 6 increase in the number of 
RLO-HMXBs at our bright luminosity cutoff, and we 
now see a non-negligible number of HMXBs at the ULX 
cutoff. At Z=0.2 Zq, we see an factor of 3 increase at our 
bright cutoff, and a factor of 30 increase at our ULX cut- 
off. At Z=0.02 Zq we see a 20% increase at our bright 
cutoff and a 40% i ncreas e at our ULX cutoff. While 
iTaam fc SandauistI ()2000D make a convincing case that 
HG-CE events should lead to mergers, the sensitive de- 
pendence of bright and especially ULX HMXB number 
on the onset of survivable CE phases demands further 
studies into this parameter. 

Another large theoretical uncertainty in the physics 
of CE phases involves their efficiency in transferring the 
potential energy from the binary orbital separation into 
ejecting the envelope of the donor star. In Figure [9l we 
show the numbe r of brigh t RLO-HMXBs as a function 
of a (defined in iWebbinlJ ([1984 ')). We note that RLO- 
HMXB number is highly conserved at all metallicities for 
all choices of a > 0.2. We find that while changing the 
value of a may determine whether specific RLO-HMXB 
progenitors survive the CE phase, the parameter space 
of RLO-HMXB progenitors is filled with "near miss" sys- 
tems (i.e. systems slightly too close to survive the CE 
phase, or systems slightly too far apart to enter RLO 
after the CE phase) such that the number of systems 
which become RLO-HMXB remains relatively constant 
with changes in the CE efficiency. In this work we do not 
consider CE efficiencies greater than unity, as no realistic 
physical mechanism has been proposed to allow for en- 
ergy loss greatly in excess of the energy available from the 
gravitational potential. We note that these Common en- 
vel ope efficiencies are in agre ement with those predicted 
bv ITaam fc SandquistI ([2000[ ). 

4.2. SG-HMXB Pathway 

From Figure[S]it is clear that SG-HMXBs occupy loca- 
tions in both observable and initial parameter space that 
are distinct from the RLO-HMXB population. While 
the systems have a luminosity function which is similar 
to RLO-HMXBs, these systems have typical orbital peri- 
ods greater than 1000 days, and occasionally longer than 
10^ days at high metalhcity. While a significant fraction 
(« 50%) of the systems have circularized due to tidal in- 
teractions, there is also a significant fraction of systems 
with moderate eccentricities. 

This observable parameter space is reasonable consid- 
ering the nature of the SG donors. Since wind accretion 
is much less efficient than RLO at moving mass onto 
a CO, only the strongest donor winds can provide the 
necessary mass loss rates to achieve bright HMXB ac- 
tivity. While highly non-spherical winds may create sig- 
nificant wind accretion with lower mass loss rates, the 
necessary spherical wind strengths are provided only by 
(super)giant systems moving though the SG phase. How- 
ever, as shown in Figure [71 (super)giant systems massive 
enough to evolve before 10 Myr have extremely large 
radii (> 1000 R©), which impose a strict lower bound on 
the orbital period of the SG-HMXB. 

The necessary orbital separations for SG-HMXBs re- 
quire a specific evolutionary chanel for their progenitors. 
Unlike RLO-HMXBs, these systems do not move through 
a CE before the SN of the primary star. A CE driven 
pathway would create much smaller orbital separations 
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Fig. 10. — Average number of SG-HMXBs with Lx > 
10^'' ergs~l over the first 20 Myr after a starburst of 10^ Mq for 
different values of the maximum fallback fraction for black hole su- 
pernovae. (Maximum fallback fraction varies from default model, 
lOx Eddington hx allowed) 

(?» 10 Rq) than the necessary separation for the SG- 
HMXB, causing the system to merge before the SG donor 
has formed. Instead, most SG-HMXB progenitors be- 
gin with a large periastron separation (> 1000 R© at 
Z:=0.02 Zq and > 5000 Rq at Z^Zq) and undergo no 
significant binary interactions before the primary SN. A 
possible exception to this rule is the case of light mass 
transfer for binaries with mass ratios close to unity. 

Since the eventual system must have a significant or- 
bital separation at the time of CO formation of the pri- 
mary, if CO formation were accompanied by an energetic 
natal kick (with a velocity > 50 km s~^), the vast ma- 
jority of SG-HMXB progenitors would be disrupted by 
the SN. For young systems, the only method for avoid- 
ing binary disruption due to natal kicks is to create the 
CO through direct collapse into a massive BH, which our 
simulations assume to impart no natal kick to the sys- 
tem. We find that this avenue produces the majority of 
our SG-HMXBs. 

In Figure [TOl we show the number of bright SG- 
HMXBs as a function of the maximum fallback frac- 
tion for BHs. Since BHs in our simulations receive kicks 
scaled down by {1-i fallback), this effectively enforces the 
assignment of a non-zero kick to all BH events. This 
shows not only that SG-HMXBs form almost exclusively 
from direct collapse BHs, but that the survival of these 
systems is radically changed if even small kicks are ap- 
plied. This restriction is most evidence at high metallic- 
ities, where the initial binary separation is largest. We 
find that adding a 20% kick to direct collapse events only 
kills about 60% of systems at Z=0.02 Z©, but destroys 
96% of systems at solar metallicity. Because of this ef- 
fect, we note that these systems should be found close to 
the star formation region, since there is no natal kick to 
impart a large spatial velocity to the binary. 

We note that there only exists a short peak (<5 Myr) 
of bright SG-HMXBs. This occurs for three reasons: (1) 
each SG-HMXB exists only for a short time, because the 
donor star moves through the SG phase very rapidly (< 



1 Myr), (2) nearly all COs are formed before 6 Myr, 
while direct collapse events are still possible. Given the 
flat initial mass ratio distribution, we expect most sec- 
ondary stars to have similar masses to the direct collapse 
primary, and thus move through the SG stage shortly 
thereafter, (3) for very small donors, the constraints on 
the radial separation from the donor star will create sys- 
tems too dim to be detected above our luminosity cutoff. 

The metallicity dependence of our SG-HMXB popu- 
lation is quite different than that of the RLO-HMXB 
population. The SG-HMXB population has many sys- 
tems which fall below our minimum luminosity cut for 
bright HMXBs. Thus, metallicity-dependent changes in 
HMXB number occur primarily due to changes in the 
luminosity of SG-HMXBs. The number of bright sys- 
tems peaks at our intermediate metallicity of Z=0.2 Zq. 
This is due to two competing effects: (1) stronger winds 
for high-metallicity donors, and (2) smaller HG radii for 
low-metallicity donors, allowing SG-HMXBs to form in 
closer orbits. 

Finally, we note one important parameter space effect 
which may greatly affect the statistics of our population. 
Here we have shown only the orbit-averaged luminosities 
of our SG-HMXB population - while we know that the 
bright and ULX HMXB population usually exhibits some 
luminosity variability. From Figure [SI we note that be- 
tween 40%-80% (for Z=Zq and Z=0.02Zo respectively) 
of our SG-HMXB population has eccentricities smaller 
than 0.05, which corresponds to a luminosity variation 
of less than 22% from periastron to apastron. At higher 
resolution, we note that most of these eccentricities are, 
in fact, almost identically zero, due primarily to tidal 
effects. However, our remaining SG-HMXB population 
has a nearly fiat eccentricity distribution, which creates 
a subpopulation of SG-HMXBs with potentially sizable 
luminosity variations (at least a factor of 50 change for 
the 10% of solar metallicity HMXB with eccentricities 
greater than 0.75). 

There are two effects which limit the impact of HMXB 
variability on our results - especially in the ULX regime. 
First, in Figure [21 we note that the XLF of our source 
population is extremely hard. Since HMXBs will only be 
near periastron for a small fraction of their lifetimes, we 
would need a large underlying population of HMXBs just 
under our ULX cut in order for variability to cause these 
systems to make substantial transient contributions to 
our ULX population. Secondly, in the same Figure, we 
note that the creation of our brightest HMXBs is con- 
trolled primarily by our allowance of super-Eddington 
accretion. Thus periastron orbits in our brightest sys- 
tems will not create significantly brighter HMXBs than 
we already see - as this luminosity will be capped by our 
luminosity limit of ten times the Eddington luminosity. 
However, we note that this variability may have an ef- 
fect at much lower luminosities, and it may create some 
variability in the number of observed SG-HMXB sources. 

5. DISCUSSION 

Our simulations yield a number of key results: (1) 
the young, extragalactic HMXB population is domi- 
nated by systems moving through two specific path- 
ways, RLO of main sequence donors onto reasonably 
sized BHs, and wind accretion from (super)giant donors 
onto massive BHs formed via direct collapse (consis- 
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tent with iVan Bever fc Vanbevereiil (|2OOO0 ): (2) while 
metaUicity does not greatly affect the characteristics of 
HMXBs moving through either pathway, it affects the 
relative strength (in number of systems) of each path- 
way, creating a total HMXB population which looks 
very different at different metallicities; (3) with mild al- 
lowances for super-Eddington accretion, the majority of 
the ULX population can be explained as an extension 
of the stellar HMXB population (consistent with both 
observations (Irwin ct al 2004: 'Pakull fc Mirionil 120021 : 
[Zampicri ct al. 2004; Swartz ct al. 2 008) and theoretical 
models jM adhusudhan ct alj 120081 : iMapelli et all 120091 : 
iZampierifc Robcrts.2009a )'): (4) as opposed to previous 
models, we find this phenomenon to be connected to the 
binary evolution of HMXB progenitors, rather than to 
the fact that low metaUicity stars produce more massive 
BHs. 

5.1. ULX Population 

Re-examining our earlier analysis, we note several ob- 
servable results in our description of the ULX population 
due to stellar HMXBs. First, in Figure [U we see that 
the 0.02 Zq environment favors ULX-HMXB formation 
by approximately a factor of 5 over solar metaUicity en- 
vironments, with a variation of approximately 2 between 
0.02 Tjq and 0.2 Z©. However, this relation is highly 
time dependent, and higher metaUicity environments are 
favored for very young star forming regions (<10 Myr). 
Since the high metaUicity population is dominated by 
SG-HMXBs, while the low metaUicity environment fa- 
vors RLO-HMXBs, in Figure[Sl we note that high metal- 
licity ULXs should have substantially longer orbital pe- 
riods than low metaUicity systems (see Figure |4]). 

There are several reasons to believe that our simulated 
ULX population robustly describes the observed popula- 
tion. In Figured (top), we note a sudden cutoff in HMXB 
luminosity at the Eddington limit. This dropoff in- 
cludes nearly 30% of our sources with Lx > lO'^^ergs^^, 
and is not observationally detected. Instead, our use 
of mild super-Eddington accretion (bottom), allows a 
much smoother decline, which is maintained up to lu- 
minosities of 10^° erg s~^, consistent w ith observations 
of nearby galaxies (jGilfanov et al.ll200^ . Secondly, using 
this method, we can replicate a large percentage of the 
ULX population, with a reasonable percentage of sources 
as bright as ix w 10^°ergs~^. 

Recent observations have found an inverse correlation 
between clu ster metaUicity and ULX fo rmation rates 
(lSoriall2d07l ). It has recently been posited (|MaDelli et al.l 
l2009f) that the reason for this trend is that low metaUicity 
stars lose less mass to stellar winds, and are thus able to 
produce significantly more massive COs than high metal- 
licity stars. On its face, this hypothesis is reasonable, 
since the primary restriction on forming ULX systems 
from the stellar population comes from the Eddington 
limit, and the Eddington limit increases linearly with 
CO mass. Thus massive BHs (« 30 Mq) would require 
only mild violations of the Eddington limited luminosity 
in order to explain the majority of the ULX population. 

While our own modeling codes agree with the corre- 
lation between low me taUicity and large BH mass (see 
iBelczvnski et al.l l2006f) , there are several problems in 
correlating this r esult to explain the U L X population, 
as done in .Madhusudhan et all (|2008D : IMapelli et al.l 



(|2009D : [Zampieri fc RobertsI (|2009aD . It is important to 
note that ULX-HMXBs form from only a very small per- 
centage of the initial binary population, and are thus 
likely to have properties quite different than the aver- 
age binary (or average HMXB for that matter). While 
it is very easy to form large BHs in low metaUicity en- 
vironments, it is actually quite difficult to find them in 
HMXBs with ULX luminosities for two reasons: (1) the 
SG-HMXB pathway is unable to create ULX HMXBs 
because the stellar wind from the low-metallicity donor 
stars are too weak to achieve these luminosities, and (2) 
it is difficult to create massive BHs through CE phases, 
since this tends to strip a large percentage of the primary 
envelope. 

5.2. NGC 1569 

A second interesting result from our models per- 
tains to the HMXB population of NGC 1569, which 
is thought to be dominated by clusters less than 



20 Mvr old ([Hunter et al.|[T98l llsrael fc de Bruvn|[T98l 
Wall en 11991 1 . Usi ng the Chandra X-Ray telescope, 
Marti n et al.l (j2002[ ) found 14 X-Ray point sources in 
NGC 1569, 1 2 of which are probable X-ray binaries. 
iKaaret et al.l (j2004[ ) have found that these X-Ray bina- 
ries are found preferentially close to the largest starburst 
cl usters in NGC 1 5 69. 

lArp fc Sandagel (jl985[ ) found two massive starburst 
clusters, aptly named Clusters A and B. Both clusters are 
thought to have ages o f less than 2 Myr. We note tha t 
Ide Marchi et al.l (|1997D . and later lOrigha et all (|200lD . 
used WFPC2 observations to conclude that Cluster A is 
actually composed of two separate clusters of slightly dif- 
ferent ages, a res ult that is not taken into account here. 
Later studies by iHunter et al.l (|20 0Q) and lAnders et al.l 
(J2004D found 48 and 169 clusters respectively. However, 
the majority of these clusters have masses several orders 
of magnitude smaller than Cluster s A or B. 

According to the estimates of lAnders et al.l (|2004[ 1. 
both Clusters A and B have ages of ~12 Myr, and 
masses of 1.64 x 10^ M© and 5.65 x 10^ Mq respectively 
- though recent studies have shown t he Cluster B m ass to 
perhaps be as much as 10% higher (iLarsen et al.|[2008i) . 
Thus, they would naively be expected to both contribute 
a significant fraction of the bright HMXB population. 
However, our results cast an interesting light on the re- 
lationship of Clust ers A a nd B to the resultant HMXB 
population. Anders et al.l ([2004) find Cluster A to have 
a metaUicity of Z=0.02 Z©, while Cluster B has a much 
higher metaUicity of Z=0.4 Z©. From Figure [1] at a lu- 
minosity cutoff of 1 X 10^^ erg s~^ (top), the Cluster A 
population is preferred by approximately a factor of 3. 
When accounting for the mass ratio of the two clusters, 
which is approximately 3 to 1, we find that the HMXB 
contribution from Cluster A is about 9 times as large as 
that of Cluster B, meaning that it will almost entirely 
dominate the HMXB population of NGC 1569. We note 
that this result would be much different, if Clusters A and 
B were instead only 8 Myr of age, and the significant SG- 
HMXB population of Cluster B was still non- negligible. 

6. CONCLUSION 

We have studied the young, extragalactic, HMXB pop- 
ulation and determined that it is composed of systems 
moving through two distinct binary pathways: mass 



11 



transfer systems main sequence donors and wind accre- 
tion systems containing (super)giant donors. Metallic- 
ity greatly affects the fraction of systems which move 
through each pathway, but generally does not change 
the characteristics of binaries moving through a specific 
pathway. From these results, we have made three ob- 
servationally testable claims: (1) the majority of the 
ULX population can be explained as an extension of the 
stellar HMXB population, (2) ULXs should be preferen- 
tially found in low mctallicity clusters - while the ULXs 
found in higher metallicity regions should be younger and 
have larger orbital periods, and (3) the starburst galaxy 
NGC 1569 is dominated by contributions from Cluster A, 
which likely provides approximately 90% of the observed 



HMXB population. Detailed modeling of both theoreti- 
cal phenomena, such as the effect of non-spherical winds 
falls outside the scope of this paper, but is currently on- 
going. New observations of starbursts in low metallicity 
environments are also underway, and will have the power 
to quantitatively test our predictions and determine the 
important role that metallicity plays in HMXB evolution. 
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